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It has been speculated that gravity could be an emergent phenomenon, with classical general 
relativity as an effective, macroscopic theory, valid only for classical systems at large temporal and 
spatial scales. As in classical continuum dynamics, the existence of underlying microscopic degrees 
of freedom may lead to macroscopic dissipative behaviors. With the hope that such dissipative 
behaviors of gravity could be revealed by carefully designed experiments in the laboratory, we 
consider a phenomenological model that adds dissipations to the gravitational field, much similar 
to frictions in solids and fluids. Constraints to such dissipative behavior can already be imposed by 
astrophysical observations and existing experiments, but mostly in lower frequencies. We propose a 
series of experiments working in higher frequency regimes, which may potentially put more stringent 
bounds on these models. 
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I. INTRODUCTION 

Gravity has not yet been unified with other fundamen¬ 
tal forces, which have been put into the Standard Model 
of particle physics. General relativity, currently the most 
successful theory for gravity, describes gravity as the clas¬ 
sical geometry of space-time, which interacts only with 
the classical energy-momentum content of matter. No 
fully consistent quantum theory of gravity has been found 
— nor is it completely clear how quantum gravity can be 
experimentally probed. 

General relativity has been systematically verified with 
increasing accuracy in the classical domain, with more 
opportunities arising from higher precision laboratory ex¬ 
periments and gravitational-wave detection [1, 2[. On the 
other hand, having achieved great success in the micro¬ 
scopic world, quantum mechanics are now going to be 
tested in systems involving macroscopic mechanical ob¬ 
jects, where gravity may start playing a role [3-14]. 

In this paper, we shall still start off as testing general 
relativity, namely, testing the speculation that general 
relativity is not a fundamental theory, but an emergent 
one that appears to take place only in macroscopic spa¬ 
tial and temporal scales. In this work, we propose an 
experimental strategy to search for evidence of emergent 
gravity. We will soon see that the tests here will be in¬ 
timately connected to tests of quantum mechanics, in 
systems where gravity plays a role. 

Emergent gravity has been motivated from at least 
two types of reasoning [15] . The hydrodynamical point 
of view, proposed by Sakharov [16, 17], argued that Ein- 
steinian gravity can naturally arise at large scales due to 
the dependence of vacuum energy (of non-gravitational 
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fields) on spacetime geometry. The thermodynamical 
point of view focuses on the role played by horizons, ini¬ 
tially from the fact that the area of an event horizon 
never decreases [18] and the creation of particles by black 
holes [19]. Later, the work of Jacobson made a more gen¬ 
eral extension to Rindler horizons of particles that can 
exist anywhere in a spacetime [20], and made the con¬ 
nection between Einstein’s equation and an equilibrium 
equation of state. Recent work further considered the 
microscopic origin of the horizon entropy — motivated 
by the holographic principle [21] and by loop quantum 
gravity [22]. 

A macroscopic dynamics emerging from an averaging 
over microscopic dynamics often has two features: (i) 
fluctuations can arise from the imperfect averaging over 
the microscopic degrees of freedom, (ii) the macroscopic 
equations of motion can be non-conservative (often time 
asymmetric, or dissipative) due to excitations of the mi¬ 
croscopic degrees of freedom that do not instantaneously 
couple back to the macroscopic dynamics. An instruc¬ 
tive example is the emergence of classical continuum dy¬ 
namics from molecular dynamics and (quantum) atomic 
physics. The macroscopic properties of solids and liq¬ 
uids, for example: (i) density, elastic moduli, compress¬ 
ibility, as well as (ii) thermal conductivity, friction coeffi¬ 
cient, viscosity all in fact emerge from the microphysics of 
molecules. Because fundamental physical laws are usu¬ 
ally considered to be time symmetric, those properties 
in group (ii), which are associated with time asymmetric 
physical processes, are direct signatures of the underlying 
molecular dynamics. In fact, physical processes associ¬ 
ated with these properties had been discovered and well 
characterized long before the underlying microphysics 
was well understood — which took place after more di¬ 
rect observations of molecular structure became feasible. 
The above analogy generates hope that even though the 
microphysics that underlies gravity may well be taking 
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place at very high energies and very small spatial and 
temporal scales, astrophysical observations and labora¬ 
tory experiments in the currently accessible regime might 
be able to reveal their existence. 

Even without specifically considering the concept of 
emergence, testing the validity of general relativity has 
long been a major research effort [1, 2]. However, most 
of these tests were performed in astrophysical or cosmo¬ 
logical settings, for example in the solar system [23, 24] 
and in relativistic binaries [25, 26]. In the context of 
laboratory tests, deviation from Newton’s inverse square 
law at small distances has been proposed as a possible 
signature of models with extra dimensions [27-29]. Ex¬ 
periments have been performed to put bounds on such 
deviations [30]. Experimental bounds on Lorentz invari¬ 
ance (of many kinds) have also evolved substantially [31- 
34] since the days of the Michelson-Morley experiment 
and now tightly constrain several alternative gravity the¬ 
ories [35, 36]. All such tests focus on low frequencies 
— since the theories being tested for mainly have conse¬ 
quences in these regimes. In this work, we shall consider 
families of theories motivated by emergence, in particular 
those that have signatures toward higher frequencies and 
smaller length scales. We will consider a variant of non¬ 
local models [37, 38], which were indeed motivated by 
space-time having “memories” in the past. In these mod¬ 
els, the gravitational field lags the motion of the source, 
even in the near zone. We will focus on tests on spatial 
and temporary scales accessible in the lab. 

From a different starting points, experiments have also 
been proposed for testing quantum mechanics involving 
macroscopic objects [3-14]. Some of the theoretical mod¬ 
ifications of quantum mechanics were motivated by grav¬ 
ity [5, 9, 39, 40], some are motivated by the determinism 
of the classical world [41, 42], some both [39, 40]. Most 
of these experiments focused on the stochastic aspect of 
the modifications, although dynamical effects were also 
speculated [5, 43, 44], 

In this work, we will recognize that our variant of non¬ 
local gravity [37, 38] and the Diosi-Penrose (DP) gravity- 
induced decoherence [39, 40] can be unified — both can 
arise from an infinite family of fields that couple with the 
energy-momentum content of space-time: non-local grav¬ 
ity as the dynamical consequence of the coupling, while 
DP decoherence as the stochastic back action. Such a 
connection has been hinted by Diosi [43, 44], but we shall 
make a more general argument. We will then explore 
possible experimental configurations that will test this 
unified model. 

This paper is organized as follows: In Section II, we 
briefly explain how dissipation might be incorporated 
theoretically, and set up a phenomenological model to 
test in the weak-field regime, which is parametrized by 
phase-lag (see the discussion in Sec. II A), similar to 

the loss-angle used to characterize material losses [45]. 
In Section IV, we briefly discuss how astrophysical pro¬ 
cesses might be used to constrain </>(w), and argue that 
there exists a new regime that is particular suitable for 


lab experiments. In Sec. V, we propose experimental 
strategies and explore their potential performance. 

II. DISSIPATION IN SPACE-TIME 

In this section, we set up a strawman theory to fit 
in our intuition about spacetime dissipation. In partic¬ 
ular, we would like to explore dissipative modifications 
to general relativity where the dynamics of spacetime is 
changed. This could be due to tracing out unknown “en¬ 
vironmental" or microscopic degrees of freedom, but here 
we focus only on a effective theory for the spacetime itself. 
In this sense, these strawman models are phenomenologi¬ 
cal, as the unknown microphysics no longer enters the dy¬ 
namics, but rather affects the parameters of these models. 
It is also possible that different micro-physical theories 
lead to similar phenomenological models at low energies 
and large spatial scales. 

In addition, we expect dissipative effect leads to a dis¬ 
tinctive experimental signature in the linearized gravity 
regime: Newton’s constant will effectively become com¬ 
plex in the frequency domain, so that for an oscillatory 
source, the phase of the Newtonian gravitational poten¬ 
tial will lag behind that of the oscillatory of the motion, 
with a finite phase difference, even in the near zone. We 
shall discuss the phenomenology of such models, illus¬ 
trating how they differ from traditional modified gravity 
models. 


A. Non-local Einstein’s Equation and its relation 
to existing modifications to GR 

Let us assume that the presence of unknown micro¬ 
physics introduces a delayed-response for the spacetime 
geometry with respect to its matter source. This delayed- 
response effectively makes the Einstein equations to be 
non-local, which may be written as 

/ +oo 

K^ a p{t-t')T aP (t',yi)dt'. (2.1) 

-OO 

This is a generalization of non-local gravity models pro¬ 
posed by Hehl and Mashhoon [37, 38]. The form of 
Eq. (2.1) is already written into a 3+1 form, which indi¬ 
cates that the Einstein tensor at one space-time event is 
determined as a weighted average over the worldlines of a 
family of fiducial observers. The existence of these special 
observers implied a preferred frame — as a consequence 
of the unspecified miscrophysics, similar to the effect of 
extra degrees of freedom in other modified gravity theo¬ 
ries — for example the Einstein-aether theory [46] , where 
the configuration of the aether field effectively defines a 
preferred frame. 1 


1 Some other approaches to nonlocal gravity include [47, 48] lion- 
commutative geometry effects. 
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The filter K lll/a p is a bi-tensor , and parametrizes pos¬ 
sible delayed-responses. It is symmetric in its first and 
second pairs of indices (i.e. K^ a/5 = A0 /tQ/3 = K, lu p a ), 
in order to respect the symmetry of G and T^ v . Since 
we have already chosen a preferred set of observers for 
the non-local integral, let us also perform a 3+1 decom¬ 
position of different components of the tensor To 

illustrate properties of such models, we choose a filter 
function of 

f 1 ) 

+ [K(t - t') - S(t - t')]5a05/30ff M V° 

( 2 . 2 ) 

so that only the 00 component of the Einstein equations 
is changed. In the vacuum case we recover the source- 
free Einstein equations G = 0, which suggests that the 
propagation of gravitational waves is unchanged in this 
modified gravity theory. Therefore in order to experimen¬ 
tally test its effect, we need to study the cases where the 
stress energy tensor of matter is nonzero. As explained 
in Sec. IIC, the matter equation of motion is modified, 
which can be derived from Eq. 2.1 and the Bianchi iden¬ 
tity. 

In addition, Eq. 2.1 shows that the Einstein tensor G^ v 
depends non-locally on the stress-energy tensor in 
other words, space-time curvature depends non-locally 
on its energy-momentum content. In order for this de¬ 
pendence to be causal, we must impose 

K(t) = 0, t < 0. (2.3) 


which leads to 

K(t) = -e~ t/T *0(t), (2.8) 

T* 

with 0 the Heaviside step function. With this filter func¬ 
tion, the space-time geometry has a (short) response time 

of T*. 

We assume that the dynamics of microscopic degrees 
of freedom only takes place at high frequencies, while in 
low frequencies accessible by astrophysical processes and 
our experiments, <j> can be simply Taylor expanded as 

0( w ) = w0'(0) + 4> r (0) + ... . (2.9) 

Since the leading order term of the above expansion is 
the same as Eq. 2.7, we can identify 

r*«0'( 0)«^. (2.10) 

U) 

which all have the physical meaning of time-lag. 

B. Linearized Einstein’s Equation and its solution 
for periodic source 

In the case where gravity is weak, we can simplify 
Eq. 2.1 by considering a perturbed flat metric = 
with \h\ <C 1. Within such limit, Eq. 2.1 can be 
linearized, while the only different component from lin¬ 
earized Einstein’s equation is (C.f. Eq. 2.1 and Eq. 2.2) 


Let us write the Fourier transformation of AT as 



dt e iult K(t) = e ^ , 


(2.4) 


where 0(cu) is expected to be small. If 0(w) is a rational 
function, condition Eq. 2.3 can be achieved by requiring 
0(cu) to have no poles on the upper half complex plane. 
Furthermore, in order to keep G and 20 v real-valued 
in the time domain, we need to impose K(t) G R, or 


[-uj 2 + V 2 ] h 00 = -16tt e^T 00 . (2.11) 

Here h is the trace-reversed spacetime metric pertur¬ 
bation, 

hfiv = 5 (2.12) 

which satisfies the gauge condition (with 4-dimesional 
covariant derivative) 


0(w) = —0*(— ui) , cuGR. (2-5) 


V' , V = 0. (2.13) 


In particular, 0(0) G R, which can be absorbed into the 
definition of the Newton’s constant, leading to 

0 ( 0 ) = 0 . ( 2 . 6 ) 

Furthermore, in order for 0 to represent a phase lag 
(instead of a variation in the magnitude of the New¬ 
ton’s constant), we shall require 0(cu) to be real-valued 
for low frequencies — although in general requiring 0(cu) 
to be real-valued for all frequencies may conflict with the 
requirement that K be causal. One prescription is to 
choose 


Note that in order to obtain Eq. (2.11) we have performed 
a 3+1 split of coordinates, transformed t into the fre¬ 
quency domain, and used x to denote spatial coordinates. 
As discussed in Sec. II A, 0(w) is a real-valued function 
characterizing phase-lag as functions of frequency, and 
7)^(0!, x) is the stress-energy tensor of the source. The 
gauge condition described in Eq. (2.13) considerably sim¬ 
plifies the form of linearized Einstein’s equation, and it 
is still compatible with dissipative-gravity modifications, 
as the left hand side of Eq. (2.1) is unmodified. 

We now solve the wave-equation Eq. (2.11) and obtain 
a retarded solution 




1 

1 — iwr* ’ 


or 0(w) « cut* , 


r* > 0 (2.7) 


^00 (/,x) = 4 


-Lpo (/, X ) i(27r/|x— x’ 

|x-x'| 


(2-14) 
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Suppose the source is oscillating at a constant frequency 
/o = u>o/2tt. In the near zone (2irfoR/c <C 1) we have 

Mf,x) ~ 4 [ d 3 x' (2.15) 

J |x-x'| 

Again, the form of Eq. (2.11) together with the require¬ 
ment that and be real-valued in the time do¬ 
main, require that </>(0) = 0. The power series expan¬ 
sion in Eq. 2.9 then allows us to switch freely between 
(f> and <f)' at low frequencies. Physically cf> carries the in¬ 
terpretations of a phase lag while <fJ is a time lag. We 
also note that this expansion does not necessarily hold 
at high frequencies, where microphysics of the underly¬ 
ing “environment" may contribute to a very different <j>. 
For concreteness, subsequent sections detailing existing 
constraints and experimental sensitivity, which are sup¬ 
posedly targeting much lower frequencies comparing to 
underlying microscopic motions, will have final results 
phrased in terms of (f>. 

To see how this model modifies near zone interactions, 
we work with the gravitational potential 1 /dg of our 
model, defined as 

floo = Voo + ^oo (2-16) 

= —1 + 2{7 DG - (2.17) 

Using this with Eq. 2.15 and Eq. 2.12, for a periodic 
source moving at frequency /o, we see that 

Udg = / dx' (2.18) 

J |x-x'| 

where p(f, x') = Too(t,x'). Note that this differs from 
the usual Newtonian potential, U, simply by a phase, 
i.e., t/ DG = Ue i ^ u °\ 

C. Bianchi Identity and equations of motion in the 
Newtonian limit 

Now let us consider the equation of motion for matter 
based on Eq. (2.1). Unlike last section where we focus 
on the metric in the wave-zone, here we are mainly inter¬ 
ested in the matter motion in the Newtonian near-zone. 
The Bianchi identity = 0 requires that 

G liV tV AY li va G av +T v va G l “ x =0, (2.19) 

where the Einstein tensor Gis to be replaced by the 
right hand side of Eq. 2.1, and the above equation be¬ 
comes the modified equation of motion for matter. Be¬ 
cause of the modification, the matter motion is generi- 
cally non-conservative. As a simple example, applying 
Eq. 2.19 and Eq. 2.1 to a point mass in the Newtonian 
limit, the equation of motion is just 

/ OO 

dt'K{t — t')T 00 {t', x) 

-OO 

1 • f°° 

~m |V(f)] 0 — m-d l h 00 / dt'K(t — t'). (2.20) 

^ .1 — 


As we anticipate that the theory recovers the Einstein 
gravity in the static limit, /jj° K{t — t')dt' should be set 
to one. As a result, the point mass equation of motion 
becomes 

a l {t) = ® , ( 2 . 21 ) 

m 

where /* is the Newtonian gravitational force acting on 
the point mass. As we can read from Eq. 2.1 and discus¬ 
sion in later sections, the metric is affected not only by 
local matter stress-energy, but also contributions of the 
past light-cone. In particular, 

^ = [ didK{t - t')^ At', x(t')), (2.22) 

m J -oo 

and here 4' is the Newtonian potential. Comparing the 
above equation with the point mass equation of motion, 
we find that 

/ OO 

dt'K(t — t')^ :i {t', x(t')). (2.23) 

-OO 

This means that within the Newtonian limit, this dis¬ 
sipative gravity model leads to an different equation of 
motion from General Relativity. 

D. Relation between non-local gravity and 
spontaneous wave-function collapse 

Let us turn to collapse models of quantum mechanics. 
The initial motivation for these models were from the 
randomness of quantum-state reduction, and its incom¬ 
patibility with classical determinism. The assumption 
of these models are that quantum states of macroscopic 
objects spontaneously become localized due to an intrin¬ 
sic collapse process [39-42] . As was later realized, these 
collapses can all be modeled in general as a continuous 
measurement on matter density [8, 10] — with spatially 
distributed measuring devices that have different mutual 
correlations. 

In a continuous quantum-measurement process, back 
action takes the form of back-action noise, but also some¬ 
times in the form of dissipation. This has lead Diosi to 
further propose that the collapse models may also cause 
the Newtonian gravitational potential to have a delay r* 
when responding to matter density changes [44]. This is 
the same as Eq. (2.9) as we consider the lowest-frequency 
contribution from non-local gravity. Diosi argued that 
current experimental data can constrain r* to around 
1 ms. However, in this paper we shall remain flexible 
about the function form of (f>{ui). 

Let us further argue that a more general can 

also arise from the quantum-measurement consideration. 
More specifically, a continuous measurement process can 
be modeled as coupling the observable we need to mea¬ 
sure with a field degree of freedom, which has an incom¬ 
ing state, and the out-going state contains information 
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about that degree of freedom. Quantum or thermal fluc¬ 
tuations in the incoming state provides the stochastic 
back-action, while information contained in the out-going 
state corresponds to dissipation. 

Now suppose we apply this to space-time geometry. 
If hoo(t,x) is being continuously monitored, we need to 
introduce an extra dimension r\, and an additional field 
n — for each x — that propagates on the t-r) plane. The 
field couples to hoo as ry approaches 0 — it therefore gains 
information about hoo, and acts back to the dynamics of 
/loo- In this way, depending on the propagation law of 
n along the extra dimension, and the detailed way it is 
couple to hoo, arbitrary shapes of the phase delay (j>(u>) 
between hoo and energy density can be constructed. 

In general, as we consider different x, coupling to dif¬ 
ferent components of the metric, as well as the above 
dependence on 77 , we can recover general non-local theo¬ 
ries of gravity. 


III. EFFECTS ON ONE AND TWO-BODY 
MOTIONS 

In this section, we shall discuss the dynamical effect of 
emergent gravity on one and two bodies. 


Writing the Fourier transform of i(f>(uj) as <E>(f) (which is 
K{t) — <$(£)), the time-domain version becomes 


Goo = 8n 


L oo 


$(t - t')T 0 o(t')dt' 


(3.4) 


Assuming $(£ — £') is a causal kernel implies <&(£ — £') = 0 
for t <t ', allowing us to extend the limits of the integral 
to ± 00 , resulting in our final expression 


Gqo = 87r 




$(/ - t')T 0 o(t')dt' 


(3.5) 


B 



Figure 1. Binary motion taking into account the dissipative 
gravity effect. In the Newtonian limit, star A feels the instan¬ 
taneous gravitational attraction from star B, as well as the 
contribution form “past" gravitational field due to the dissi¬ 
pative gravity effect, according to Eq. (3.5). 


A. Effects on Single Body Motion 


Consider a spherical, homogeneous object (with mass 
M, radius R) moving with non-relativistic velocity v rel¬ 
ative to the preferred frame where Eq. (2.1) holds. Ac¬ 
cording to Eq. (2.23), there is a force generated by the 
“past” gravitational field of the same object. This self- 
gravitational force is dissipative for the object’s motion, 
as it generates acceleration anti-parallel to the direction 
of motion: 

CM r°° 

a = -J^ v j o (3.1) 

where G has been restored in this equation. 

If we assume that the kernel function is described by 
Eq. (2.8), the acceleration evaluates to 


GM GM , 

a = —^-r*v = ——(0)v . 


R 3 


R 3 


(3.2) 


The phase-lagging feature of the nonlocal Einstein equa¬ 
tion inevitably leads to a dissipative self-interacting grav¬ 
itational field. 


B. Effects on Two-Body Motion 


In order to calculate effects on two-body interactions, 
we take a step back and look at the “full” theory. Con¬ 
sider 


Goo(/,x)= 87 re^ 2 ^T 0 o(/,x) 

«87r[l + *0(27r/)]T oo (/,x). (3.3) 


To estimate the effect in binary motion, we consider 
two equal mass stars A and B in orbit around each other 
as shown in Fig. (1). In the Newtonian limit, star A 
and B feel the instantaneous attraction force from each 
other, which are orthogonal to their velocities. If we fur¬ 
ther take into account the dissipative gravity effect, as 
shown by Eq. (3.5), the external gravitational field act¬ 
ing on star A (and vice versa for star B) contains a piece 
that originates from the past of star B. These “past" 
gravitational fields exert a tangential force on star A’s 
orbital motion, which in turn introduces additional or¬ 
bital energy loss/gain (depending on the sign of 
with respect to the gravitational radiation. For example, 
if 4>(uj) > 0 , although the gravitational field lags behind 
the source motion, the tangential force is along the star’s 
direction of motion and and the binary motion actually 
gains energy according the positive phase lag. In this 
case we see that the rate of energy gain is given by 

M 2 

Emg ~ q— / $(t - t 1 ) sm[2nf(t - t')\dt' 

^ J —OO 

~ 9(M/) 4 / 3 0(2tt/), (3.6) 

where q = m\m^l {rn\-\-rn- 2 ) 2 is the symmetric mass ratio, 
the post-Newtonian parameter v 2 = (2t:M f) 2 / 3 , and the 
subscript MG stands for mutual gravity. As we require 
that the standard Newtonian result to be recovered in the 
static limit, this constrains the DC phase-lag such that 
<j>{f = 0) = 0. Therefore if the frequency is sufficiently 
low, on can Taylor-expand the phase and obtain 

Emg ~ qM~ 1 (Mf) 7 / 3 (f>' 0 . 


(3.7) 
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Experiment 

Upper limit on 
</>(2 tt/o) 

fo (Hz) 

Reference 

PSR J0737-3039 

7.9 x 1(T 36 

1.1 x 10~ 4 

[49] 

PSR B1534+12 

4.3 x 10“ 38 

2.8 x 10~ 5 

[50] 

PSR B1913+16 

2.4 x 10~ 39 

3.6 x 10~ 5 

[51] 

PSR B2127+11C 

1.6 x 10 -37 

3.5 x 10~ 5 

[52] 

Earth Orbital Motion 

3 x KT 17 

3.2 x 10~ 8 

[53] 

Shapiro delay 

0.014 

3.2 x 10~ 8 

[54-56] 

Torsion Pendulum 

0.1 + 

0.084 

[57] 

Quadrupole Antenna 

o.F 

60.5 

[58] 

Cantilever 

1 x 10 -3 

324.1 

[59] 

Rotor/Cantilever 

0.1 

353 

[60] 


Table I. Constraints on <j> derived from existing measurements 
at different frequencies. Quantities with a * superscript are 
estimated by assuming a 0.1 rad uncertainty in the phase at 
the specified frequency. The phase lag in the Shapiro delay 
measurement is obtained by converting the constraint on the 
time-lag estimate on the same system. 


and, denoting P b as the period of the binary motion, we 
have 


P b 2 E ' 

Therefore the accuracy of period measurement constrains 
the magnitude of possible dissipative gravity effect. We 
can then use the fact that <j>' 0 = (f>(2n f 0 ) / (2n f 0 ) to phrase 
our results in terms of the phase angle. The compact¬ 
ness of neutron stars M \^/can be determined from 
different models of nuclear equation of state, and here 
for simplicity the neutron star radius R is estimated as 
13 km. Constraints from a representative sample of bi¬ 
nary pulsars appear in Table I. We can see that pulser 
systems give really stringent upper bound on <f> at low 
frequencies (recall the 1ms constraint on r* in [44]), as 
nuclear densities are much higher than normal matter. 
On the other hand, it remains interesting to constrain (f> 
at higher frequencies, using table-top experiments. 


However, this is not the leading order dissipative grav¬ 
ity effect, which comes from the self-dissipative dragging 
force as shown in the previous section. In the binary 
system, such a damping mechanism dissipates energy at 
rate 


Esg 


2 

m\ 


ml 


+ ~F* V 2 $ 


m 


Rl 


-qMvtfo 



(3.8) 

which is about a factor of r 3 /R 3 larger than Emg- 
To help further understand this result, recall that 
quadrupole radiation of gravitational waves has power 


Ei ~ —g 2 (M/) 10/3 , (3.9) 

7r 


which is one and a half post-Newtonian order higher than 
£mg, and consequently even smaller comparing to E§q! 
A direct consequence of this observation is that binary 
pulsars should provide the best observational constraints 
on the magnitude of (f> at low frequencies, because of the 
high compactness of neutron stars. 


IV. EXISTING OBSERVATIONAL AND 
EXPERIMENTAL CONSTRAINTS 


A. Binary pulsars 


B. Solar system test 


Solar system test of gravity theories generally involves 
measuring perihelion precession, spin precession mea¬ 
surement (such as “Gravity Probe B") and test of the 
“Weak Equivalence Principal". None of the above exper¬ 
iments seems to constrain a friction-type force well [61] . 
We can nevertheless expect the period change after one 
orbital period of earth: 


AT _ 3 AE _ Q GM e T 
~ 2~E~ ~ 3 R 3 e v ’ 


(4.3) 


where M E ,R E are the earth mass and average radius 
respectively. Plugging in the numbers, we find 

AT 

(f>' 0 ~ 2 x 10 2 — s ~ 1.4 x l(T 10 s, (4.4) 


where we have used the orbital period change of earth 
due to tidal friction ~ 2.3ms/century as an estimate [53]. 
Obviously this is a much looser constraint than the binary 
pulsar tests, because normal matter density is much lower 
than the nuclear matter density. As the earth orbital 
period is longer than the periods of binary pulsars, this 
test belongs to the low-frequency regime which is ruled 
out by binary pulsar tests. 


According to our previous analysis, if we turn on the 
dissipative gravity effect, the change of fractional energy 
emission rate for a binary pulsar system is 

AE _ E sg _ ( 2tt) 2 / 3 tt(I)' 0 M / 7711 772-2 \ 

= = 32g(M/)3 [r 3 + R 3 ) 

(27t) 2 / 3 7T(/>oM 2 

~ 32g(M/) 3 i? 3 ' 1 ' 1 


C. Shapiro Delay 

Shapiro delay is the time delay of light due to the lo¬ 
cal gravitational field of a nearby massive body [56, 62]. 
For a static source, our model predicts the same static 
gravitational field as general relativity, so there is no ob¬ 
servable Shapiro delay effect. However, as we have shown 
above, this is not the case for a moving source. 
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In the solar system, the time delay of light for a mov¬ 
ing source has been measured in 2002 [54, 55], as Jupiter 
moved by the line of sight to quasar J0842+1835. The 
timing sequence of the delay signal was applied to con¬ 
strain the speed of gravity, which was measured to be 
consistent with the speed of light, with a relative uncer¬ 
tainty up to 20%. We notice that the phase-lag due to 
dissipative gravity, can be effectively translated to a shift 
in the “speed of gravity" in this case. More specifically, if 
the closest distance between the trajectory of the light is 
d and the velocity component of Jupiter moving towards 
the light trajectory is i>d, then the time lag is roughly 
bounded by 

^ < 20% , or <j>' < . (4.5) 

d 5 v d 

According to the parameters presented in [54, 55] , Vd ~ 
14 km/s and d ~ 10 6 km, and consequently we obtain 
<j>' < 7x 10 4 s. Because the moving-body correction to the 
Shapiro delay is already a second order effect, in terms of 
Vd/c, this measurement a gives much weaker constraint 
on (j)' than the test based on planetary motion. 

D. Near field experiments 

Although astrophysical tests provide the tightest con¬ 
straints on our model, laboratory-scale experiments de¬ 
signed to test the nature of near-field gravity provide con¬ 
straints at much higher frequencies. Examples include 
torsion pendulua [57] and the quadruple antenna [58], 
which look for deviations in the amplitude and do not 
report phase information. In these cases, we assume the 
experiment has observed no phase lag greater than .1 ra¬ 
dians. Other experiments involving cantilevers [59] and 
improvements thereof [60] explicitly measure the phase 
angle. Table I summarizes these constraints. 

The self-damping effect for moving bodies described in 
Sec. Ill A has the consequence that there is an intrinsic 
damping force experienced by any harmonic oscillator in¬ 
volving moving masses. From Eq. (3.2), and assuming a 
harmonic oscillator with natural angular frequency w 0 , 
the mechanical quality factor is limited by 

Q~ 1 > — T*(w 0 ) = ^(w„), (4.6) 

Wo wg 

where p = M/R 3 , and we have allowed the gravitational 
delay, </>(w), to be frequency dependent. Here we see 
that the damping effect is greater for low frequency os¬ 
cillators with dense masses. Depending on the value of 
(j), this damping could exceed other mechanical dissipa¬ 
tion effects. If this is not the case, the mechanical quality 
factor will only be slighly modified by gravitational self¬ 
damping and estimating the non-gravitational mechani¬ 
cal damping could lead to large systematic uncertainty 
in the estimation of the gravitational phase lag. For the 
oscillators described in Sec. V, we find that the contraint 



Figure 2. Plan view of a schematic diagram of the exper¬ 
iment. In this example, a rotor with 6 spokes and masses 
spins at an angular frequency to near a rigid Fabry-Perot cav¬ 
ity with an angular eigenfrequency 6 to. The laser beam passes 
underneath the masses. 

placed by direct measurement of the phase lag to be more 
stringent than the constraint set by measuring the qual¬ 
ity factor. 


E. Gravitational wave experiments 

It is clear from Sec. IIIB that our model has implica¬ 
tions for gravitational wave detection, effecting the coa¬ 
lescence of compact binaries. In general, the phase evolu¬ 
tion of the waveforms will be modified and a careful study 
of the full implications of Eq. (2.1) is required. These is¬ 
sues will be explored in future studies. However, for the 
model under consideration here, binary pulsars should 
provide a stronger constraint for low frequency motions. 


V. EXPERIMENTAL PROBE OF 
GRAVITATIONAL DISSIPATION 

A. High Frequency Eotvos Experiment 

Here we discuss one possible experimental approach for 
testing the existence of dissipative gravity, this approach 
roughly models the paradigm of the classic experiments 
of Baron von Eotvos [63, 64] while incorporating the ex¬ 
perience gained over the past several decades of short 
scale gravity experiments [30, 57, 59, 65, 66] The experi¬ 
ment consiststs of an active mechanical Attractor and a 
high sensitivity gravitational Responder. They are con¬ 
structed such that their oscillations cause tidal gravita¬ 
tional interactions. By measuring the response of the Re- 
spnder, we can look for the possible phase lag predicted 
by a dissipative gravity model. 

The outline of this Section goes as follows: in Sec¬ 
tion VA1, we estimate the magnitude of the gravita¬ 
tional acceleration induced by our periodic attractor sys¬ 
tem. in Section V A 2, we estimate the sensitivity limit 
of a mechanical responder which is limied by Brownian 










thermal noise, in Section VB, we calculate the minimum 
detectable phase lag using the system we prescribe. 


1. Tidal gravitational field of the periodic attractor 

Our aim is to probe the frequency dependence of grav¬ 
itational forces and how they may differ, in the non- 
relativistic limit, from Newtonian gravity. To this end 
we aim to use a system of moving masses that create 
a periodic tidal gravitational attraction and will use a 
sensitive mechanical system to measure the mechanical 
response to the gravitational force. An example of such 
a system is shown in Figure 2. 

The scaling of the tidal gravitational acceleration of 
the periodic attractor system is approximated by 

5 tidai = GpL, (5.1) 

where G is Newton’s constant, p is the density of the 
material used in the Attractor, and L is a characteristc 
length scale of the attractor. 


2. Acceleration sensitivity of the responder 

To maximize the signal amplitude, we propose to use a 
mechanical system with a high quality factor resonance as 
our probe of the gravitational force. The motion of the 
mechanical system will be sensed using interferometric 
metrology to maximize readout sensitivity. With such 
a narrow mechanical resonance, and in addition we will 
integrate for long times to accumulate larger signal to 
noise ratio, for a given resonator system the measurement 
is essentially single frequency. 

A detailed noise analysis will be presented in Section 
VB, and it will show that in most cases the dominant 
noise source is Brownian thermal fluctuations of the me¬ 
chanical resonator. 

Near the mechanical resonance, the force power spec¬ 
tral density of thermal fluctuations seen by the mechan¬ 
ical responder is given by 

S$=4mk B T U £, (5.2) 

'°l 

where ks is Boltzmann’s constant, m is the mass of the 
responder, T is the temperature, uiq is the resonant fre¬ 
quency, and Q is the mechanical quality factor. 

Because we want to measure the sensitivity to gravita¬ 
tional acceleration, in units of acceleration power spectral 
density, the thermal noise is 

^ h = 4fc * r Sv ( 5 - 3 ) 


3. Sensitivity to phase lag 

The behavior of the mechanical resonance forces the 
Newtonian signal component to be in the orthogonal 


Parameter 

Cantilever 

Diluted 

Cantilever 

Rigid 

Cavity 

Pendula 

Responder 

Mass 

1 m g 

1 mg 

300 g 

1kg 

Responder Q 

10 6 

10 10 

10 8 

10 s 

Lowest 

Eigenfrequency 

200 Hz 

40 kHz 

10 kHz 

2 Hz 

Temperature 

100 mK 

100 mK 

IK 

1201< 

Stored Power 

10 mW 

60 W 

lkW 

10W 


Table II. Physical parameters of proposed gravitational re¬ 
sponders. 


quadrature with respect to the signal of the periodic at¬ 
tractor. And, ignoring systematic effects which will be 
discussed in Section VE, any in-phase signal will be due 
to the dissipative gravity effect. The magnitude of the 
signal in acceleration units is 

as = (5.4) 

and this is being compared to the RMS thermal acceler¬ 
ation of the oscillator due to thermal noise, 



(5.5) 


where r is the integration time, and we assume that S g is 
approximately constant near the mechanical resonance. 
Therefore the minimum detectable dissipation phase an¬ 
gle detectable with an SNR of unity after a time r is 


H^o) 


1 UksTuo 

GpL y mQr 


(5.6) 


Given the calculation shown in Figure 3, we can 
estimate that <?tidai ~ lnm/s 2 , and assuming ojq = 
2tt x 2kHz, Q = 10 s , T = IK, m = 0.3kg, and 
r = 24 hours, the minimim detectable phase angle is 
(f> = 3.7 x 10”' radians. 


B. Experimental Designs and Sensitivity Limits 

The example given in Section V A 1 is illustrative in 
showing the general order of magnitude of an expected 
signal. For a more detailed analysis including a noise 
budget estimate, it is necessary to define a specific res¬ 
onator cavity geometry and readout scheme. 

Here we describe a few potential Responder configu¬ 
rations that are common in the gravitational-wave and 
opto-mechanics communities. All take the form of a high 
finesse Fabry Perot cavity but vary in the configuration 
and geometry of the cavity mirrors. 

Schematic diagrams of the configurations are given in 
Figure 4. Here we will briefly discuss the features of each 
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2 4 6 

Time [ms] 



4 6 8 

Time [ms] 


Figure 3. Differential acceleration (peak-peak) between the two ends of a Fabry-Perot cavity with a multi-mass periodic 
Attractor (cf. Fig. 2). As the number of masses is increased, the signal frequency increases, although the peak strain is 
correspondingly reduced. In all 6 plots, m = 0.01 kg, R = 0.095 m, and d\ = 0.12 m. 



(a) Cantilever Cavity (b) Cantilever Cavity with (c) Massive Rigid Cavity (d) Massive Pendula 

Optical Spring Dilution 


Figure 4. Schematic diagrams of four possible gravitational responders. 


configuration and provide a noise budget estimate assum¬ 
ing feasible experimental parameters. These parameters 
for each configuration are given in Table II. 

A pendulum suspended test mass configuration would 
be similar to the arm cavities used in large scale 


gravitational-wave observatories such as LIGO, however 
with reduced physical dimensions. 

In this configuration, the sensing cavity is similar to the 
rigid kind used for stabilization of lasers. A solid piece 
of single crystal silicon with optically contacted mirrors 
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Figure 5. Shown here are the detailed noise budgets for the four proposed cavity configurations. The horizontal axis denotes 
positive sideband frequencies separated from the cavity mechanical resonance. In most cases the total noise is dominated by 
thermal noise of the mechanical oscillator. In the case of the suspended pendulum, the low frequency noise is dominated by 
seismic vibrations. 


is used to measure the gravitational accelerations. 

In the past few decades, microcantilever force sensors 
have become widely deployed as high sensitivity force 
and acceleration sensors [67], especially in the biological 
agent detection fields and radiological. Microcantilevers 
have also been used for sensitive fundamental physics 
experiments [68], such as Casimir force experiments [69] 
and the searches for fifth force and extra dimensions by 
looking for short range deviations in the Newtonian 1 /r 2 
law. They are recently also candidates to measure quan¬ 
tum backaction noise [70, 71]. 


C. Components of the Noise Model 

Our noise model consists of several noise terms which 
limit the measurement of the gravitational coupling of 
the periodic attractor to the responder. Each noise term 
is explained below. 

To characterize the Responder, in particular, we are 
interested in the sensitivity to tidal gravitational accel¬ 
erations. To calculate the acceleration sensitivity, we 
first compute the displacement noise contribution from 


all terms, then this is multiplied by |x _1 /rn| to produce 
the acceleration noise. As described in Section V A 3, the 
sensitivity is maximized at the mechanical resonance, and 
the resonance is very narrow, thus we make logarithmic 
plots in Figure 5, where the origin of the horizontal axes 
are at the mechanical resonance frequencies. 


1. Force Noises 

a. Seismic Vibrations of the laboratory due to ex¬ 
ternal seismic fluctuations and nearby vibrating machin¬ 
ery will produce motion of the support point of the Re¬ 
sponder cavity. This will be mostly rejected below the 
Responder’s mechanical resonance frequency. Near the 
resonance the rejection will be imperfect due to imper¬ 
fect matching of the pendular resonant frequencies (we 
assume a mismatch of ~ 0.1%) of the two ends of the 
cavity. 

b. Newtonian Gravitational Noise Fluctuations in 
the Newtonian gravitational potential due to local 
sources (e.g. seismic, anthropogenic, air pressure, etc.) 
can limit the acceleration sensitivity due to their spatial 
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Figure 6. Comparison of Experimental limits and Theoretical Predictions. The arrows represent existing upper limits 
on the phase lag angle (cf. Table I). The shaded blue region is excluded due to pulsar observations. To give a sense of scale, the 
shaded regions at the bottom of the plot give estimates for <f> assuming a dissipation of unity at the Electro-Weak and QCD 
length scales. The solid lines are expected upper limits from the experiments described in Sec. VB, where each point in the 
line represents an oscillator resonant at that frequency. The thick purple line (from 1-30 Hz) represents the massive pendulum 
cavity. The thick orange line (from 30-3000 Hz) represents the small cantilever cavity. The thick grey line (from 1-100 kHz) 
represents the rigid spacer cavity. 


gradients across the ends of the Responder cavity [72, 73]. 
As most of these fluctuations are sourced by mechanical 
vibrations (and have a rather red power spectrum), they 
are most serious for the massive suspended cavity and 
least serious for the rigid, monolithic cavity (cf. Fig. 4). 

c. Thermal Noise Several different kinds of thermal 
noise are relevant for these experiments. Internal fric¬ 
tion of the suspension support elements, cantilever flex 
joints, and rigid cavity spacer result in Brownian forces 
which move the mirrors. The same kind of internal fric¬ 
tion in the mirror coatings also produces Brownian [45] 
strain fluctuations of the coating phase. Thermodynamic 
temperature fluctuations [74, 75] can lead to apparent 
optical cavity length fluctuations through a finite ther¬ 
mal expansion coefficient of the mirror substrate (Sub¬ 
strate Thermo-Elastic) and coatings, as well as through 
the temperature dependence of the index of refraction of 
the coating materials [76, 77]. We refer to the coherent 
expansion and refractive fluctuations in the coating as 
‘Coating Thermo-Optic’. 

d. Quantum Vacuum Fluctuations of the vacuum 
electro-magnetic field limit a variety of high precision 
measurements [78, 79]. In the particular case of a 


Fabry-Perot interferometer, the fluctuations in the phase 
quadrature limit the resolution of the optical phase mea¬ 
surement and is usually referred to as ’shot noise’. The 
fluctuations in the amplitude quadrature produce radia¬ 
tion pressure ‘shot noise’ in the circulating field within 
the Fabry-Perot cavity; this fluctuating pressure moves 
the cavity mirrors. This Quantum Vacuum noise limits 
these experiments in different frequency bands. If the 
thermal noise can be reduced, then squeezed light tech¬ 
niques can be employed to further constrain the space- 
time dissipation. 


D. Limits on complex phase angle for oscillators of 
different frequencies 

The sensitivity to tidal accelerations given from our 
noise model can be combined with the acceleration in¬ 
duced by the periodic attractor to give a final mea¬ 
surement limit of the complex gravitational phase angle. 
Generalizing Equation 5.5 for a frequency varying S g , the 
measurable limit on cf) is 
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<P\im 


1 / r O + V^avg dw 

f/att y aoj 0 —l/2T avg 


(5.7) 


where Sg is the power spectral density of the accelera¬ 
tion noise of the receiver, wo is the receiver resonance 
frequency, and g att is the differential acceleration caused 
by the periodic attractor. We estimate an acceleration 
on the order of lnm/s 1 2 for g at t, as illustrated in Fig¬ 
ure 3. The limits for the various types of oscillator, for 
various possible resonant frequencies, after 24 hours of 
integration, are shown in Figure 6. 


the scatter signal we will apply a few Lambertian 
marks to the Attractor at irregular intervals. 

4. In order to make a high precision estimate of <f >, one 
must know the mechanical admittance precisely. 
The phase angle of the admittance will be mea¬ 
sured using non-gravitational forces (e.g. with elec¬ 
trostatic force actuators and radiation pressure on 
the laser beam). The difference between the elec¬ 
tromagnetic force admittance and the gravitational 
admittance will ultimately limit the sensitivity of 
this search. 


E. Elimination of Systematic Effects 

The sensitivity shown above in Figure 6 includes only 
the limits due to Gaussian random noise. In reality, ex¬ 
periments where the signal is so weak are likely to be lim¬ 
ited by systematic effects. In these searches for a gravita¬ 
tional anomaly, we should be especially concerned about 
an anomalous signal showing up due to imperfections in 
the experimental apparatus, rather than new physics. 

There are several systematic effects including time de¬ 
lay, mechanical response limits, mis-estimation of source 
oscillator phase, clock / timing errors in the digitizer, 
timing error in the interferometric measurement of the 
source motion, etc. 

We will estimate and subtract these effects in the fol¬ 
lowing ways: 

1. Electro-magnetic pickup of the electrical driving 
signal in the readout electronics at /drive will be 
rejected since it is not at the gravitational pertur¬ 
bation frequency (/signal)- High frequency harmon¬ 
ics of the drive signal will be attenuated by usual 
grounding and shielding practices [80] . 

2. Patch charges on the Attractor or Responder can 
produce a fluctuating electrostatic force at /motion- 
Both the Attractor and Responder can be coated 
with a slightly conductive coating to minimize 
charge gradients and the Attractor will be enclosed 
within a Faraday cage to ground these fields. 

3. Light scattered from the cavity mirrors can 
backscatter from the Attractor back into the cav¬ 
ity mode to produce a signal at /drive- To mitigate 
this, we will use super polished optics for low scat¬ 
ter and blacken the nearby surfaces. To calibrate 


VI. CONCLUSIONS 

The confrontation between quantum mechanics and 
general relativity points to a modification of one or both 
of the theories. Previously, two communities have been 
tackling these problems: modification of gravity, and pos¬ 
sible modifications of quantum mechanics. In this work, 
we intend to bring up the connection between these two 
different approaches. In order to do so, we have written 
down model theories that are both be viewed as modifica¬ 
tions to GR and to QM. Experimentally, these can plau¬ 
sibly be probed by precision measurement performed in 
the laboratory scale. While it is certain that the proposed 
experimental approaches are extremely challenging, they 
should be able to, at the very least, place interesting up¬ 
per limits on these types of alternative theories of gravity 
and point us towards new physics. 
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